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ABSTRACT 

We have measured Lick index equivalent widths to derive luminosity weighted stellar ages and 
metallicities for thin and thick disk dominated regions of 9 edge-on disk galaxies with the ARC 3.5 
meter telescope at Apache Point Observatory. In all cases, the thick disks are confirmed to be old 
stellar populations, with typical ages between 4 and 10 Gyr. The thin disks are uniformly younger 
than the thick disks, and show strong radial age gradients, with the outer regions of the disks being 
younger than 1 Gyr. We do not detect any significant metallicity differences or a-element enhancement 
in the thick disk stars compared to the thin disk, due to the insensitivity of the Lick indices to these 
differences at low metallicity. We compare these results to thick disks measured in other systems and 
to predictions from thick disk formation models. 

Subject headings: galaxies: stellar content — galaxies: spiral — galaxies: structure 



1. INTRODUCTION 

Old stellar populations preserve a fossil record of a 
galaxy's early formation and evolution. In the Milky 
Way, the stellar thick disk and halo have been rec- 
ognized as the oldest stellar components and have 
been studied extensively. Stars in the MW thick disk 
are old (~8-12 Gyr), and are me tal poor when com- 
pared to local thin disk stars (.Reid &: Maiewskil Il993t 
IChiba &: Beersll2000( ). Their chemical composition shows 
that they are enhanced in a-elements compared to 
thin disk s tars, suggesting a rapid formation tiinescale 
(<1 Gyr) (Reddv et al.' '2006^; 'Brewer & Cameyl l200l 
iBensbv et al. 2003, 2005; Prochaska ct al. 200°^ 

Because stellar thick disks and halos are intrin- 
sically low surface brightness features, they have 
been observed in detail in only a handful of galax- 
ies. Thick disks have been photometrically detected 
as an excess flux at lar ge galactic latitudes across a 
range of Hubble types (iDalcanton fc BernsteinI [20021: 
Bursteinl[l97l iTsikou dil 1X9791 : Ide Griis k van der Kruitl 
19961: Ide Griis fc Peletie^ [l997l: fPohlen et all 
200^ 'van der Kruit 1984; 'Shaw fc Gilmore Il989t 
van Dokkum ct al. 1994; Morrison_ct al. 1997; W u et all 
20021: lAbe et all 119991 : iNeeser'et all 120021 ) . as well as 
detected in star counts using r esolve d stellar pop- 



ulations from HST (iSeth et all [20051: iMouldl 



Tikhonov et al. 2005; 'Tikhonov fc Galazutdinoval 



2005 



2005 



Xikhonov fc Galazutdinova 2008). 

Three general classes of formation mechanisms have 
been put forward to explain thick disks. In the first, 
a previously thin disk is kinematically heated. In this 
scenario, stars form in a thin disk and increase their ve- 
locity dispersion with time. This vertical heating can 
be ra p id, due to interaction s and mergers (iQuinn et al.l 
19931: [Walker et al.l [l99fit iVelazouez fc Whitd 119991: 



Chen et al. l l200lURobin" et al. 1 119961: iVillalobos fc Helmil 



2001 " or gradual due to scattering off giant molecular 
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clouds, spiral arms, and/or dark matter substructure 
(IVillumse nlll985l: iHan ninen fc Flvnnll20 02l: [Benson et al.l 
,2004; Hav ashi fc Ch iba 2006; Kazantzi dS et al.l [20071 1 
In the second formation scenario, stars "form thick" 
with star f ormation occurrin g above the midplane of 
the galaxy ([Brook et al.[[200^ or form with lar ge initial 
veloc ity dispersions in massive stellar clusters ([Kroupa[ 
[20021) . In the final class of models, thick disk stars 
are directly accreted from satellite galaxies. Numeri- 
cal simulations have shown that stars in disrupted satel- 



([Abadi et al.[ [20031: 


Martin et all [2004[: [Bekki fc Chiba 


2001 


; Gilmore et al. 


20021 [Navarro et all [20041 [Statlei 


1988 


; [Read et al.[ [20081). While these models were orig- 



inally developed to explain the origin of the MW thick 
disk, they should work equally well for thick disks in 
other massive galaxies. 

On the other hand, some of these mechanisms are 
likely to be less effective in lower mass galaxies which 
have lower density disks, little or no spiral structure, and 
fewer satellites hosting stars. Lower mass galaxies also 
have different formation times, environments, and grav- 
itational potentials that may also lead to variation in 
the mass, age, and metallicity of thick disks of low mass 
galaxies. 

To test these formation models, detailed comparisons 
of thin and thick disk properties are required across a 
range of galaxy masses. In particular, the relative ages 
and chemical enrichment patterns of the thin and thick 
disks should differ among these formation models. If 
the thick disk is just a gradually kinematically heated 
thin disk, there should be a smooth age and enrichment 
gradient between the two. In contrast, if the thick disk is 
formed from accreted stars wc should expect the ages and 
metallicities of the thin and thick disks to be only weakly 
correlated. We may also expect to see variations with 
the mass of the galaxies, with less massive galaxies being 
more susceptible to external heating and more massive 
galaxies being better able to tidally disrupt satellites. 

Measuring the ages and metallicities of thick disks out- 
side the MW has proved to be challenging. Ages and 
metallicities can be derived from isochrone fits to thin 
and thick disk stars resolved with HST, provided that 
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the host galaxies are suiEciently close and oriented edge- 
on to the line of sight. The systems studied this way show 
older populations at large scale heights but little verti- 
cal metallicity gradient, at least for the low mass galax- 
ies (Vc < 100 km s~^) which d ominate these samples 
(|Seth et al.ll2005l : iTikhonov et al...2005: .Mould. 200 5). 

For the systems that are farther away, only broadband 
colors have been used to estimate the ages and metal- 
licities of the thick disks. When thick disks are pho- 
tometrically detected, they typically have very red col- 
ors (B - R 1.3 - 1 51 (lYoachim fc Dalcanto^ 12006 



vations probed to ~1 scale length. However, all of these 
galaxies were fairly face-on, preventing thick and thin 
disk components from being separated. 

In this paper, we target regions of edge-on galaxies that 
are dominated by either the stellar thick disk and thin 
components. We then compare ages and metallicities de- 
rived from Lick indices both between the two components 
and from galaxy-to-galaxy. 

2. OBSERVATIONS 

The original sample of 49 galaxies was selected from 



iDalcanton fc BernsteinI l2002t Ivan DokkTim et al.Hl99l ^^"^ Flat Galaxy Catalog ( Karadrentscv et al. 1993) and 



suggestive of an old population. However, stellar param- 
eters are notoriously difhcult to derive from broadband 
colors due to the age-metallicity degeneracy in the op- 
tical colors, and a lack of IR-colors at the low surface 
brightnesses of the thick disk region. 
Some progress has been made in measuring the metal- 



sion hnes (e.g. 


Tremonti et al.ll2004l: iZaritskv et al.lll994l: 


Ivan Zee et alJ 


1998f). and in measuring metallicity gra- 



While the emission lines studies can constrain the to- 
tal chemical enrichment a galaxy has experienced, they 
tell us nothing about the underlying stellar populations 
in the galaxies. Emission lines are also unable to pro- 
vide constraints on the properties of extra-planar stellar 
populations, which are likely to be dominated by old, 
dormant stellar populations. 

To better measure the ages and metallicities of thick 
disks in a larger sample of galaxies, we turn to the inte- 
grated spectrum of these galaxies and use the Lick/IDS 
absorption line system to derive luminosity-weighted 
stellar population properties. The Lick indices were 
originally dev eloped for studyi ng older stellar popula- 
tions (Burstci n et al.l [l984: Fab er et al.lfT985) . and have 
been used extensively in an alyzing elliptical g alaxies 
and globular clusters (e.g. 



ISanchez-Blazquez et al.l |2007| V^ 



Trager et all Il998l l2000al : 
' iWorthev et all (IT991 ) 



showed that using a combination of age sensitive (Balmer 
lines) and metallicity sensitive (Mg h and Fe) indices, one 
can lift the age-metallicity degeneracy for a stellar pop- 
ulation. Stellar spectral libraries have now been used to 
create SSP models over a large ra n ge of metallicity and 
age combinations (IWorthevI 11994 IVazdekis fc Arimotol 
11999), incl uding models with variable a-element en- 
hancement ([Thomas et al.l[2"003[ ). The tools also now ex- 
ist to calcula te expected Lick indices fo r composite stellar 
populations (|Bruzual fc Charlo^l2003f ). 

Despite the development of stellar synthesis codes that 
can be extended to younger stellar populations, rela- 
tively few studies have attempted to observe Lick in- 
dices in disk galaxies. Studies using tunable filters 
have be en able to detect Mg and Fe index gra dients 
in disks (jBeauchamp fc Hardvlri997l : iRvder et al.ll2005f ). 
but have not been combined with measurements of age- 
sensitive indices. Studies of disk systems have tended 
to focus on the high surfac e brightness bulge compo- 
nents (iMoorthy fc Holtzma nI l2006t iPeletier et all 120071 : 
iPrugniel et all I2OOI, : .Proctor et al.l l2000f) . a^^fail to 
reach very far into the disks. In the m ost extensive stud ; 
observing Lick indices in disk galaxies. lMacArthurl (j200i 
observed Lick indices in 8 galaxy disks, including several 
late- type galaxies out to ^1 scale length. These obser- 



imaged in B, R, and K s ()Dalcanto n fc Bernstein 2000). 
IDalcanton fc BernsteinI (|2002| ) used this imaging to 
demonstrate the ubiquity of thick disks around late-type 
galaxies. We have since used two-dimensional decompo- 
sitions of the galaxy images to measure structural param- 
eters for the thick and thin disks (Yoa chim fc Dalc antonI 
l2006f ). We have also measured kinematic proper- 
ties of the thick and t hin disks using GMOS on 
the Gemini telescopes ([Yoachim fc DalcantonI 120051 : 
lYoachim fc Dalcantonl l2008'l 

For this study, we selected a subset of galaxies where 
the photometric decompositions suggested we would be 
able to isolate thin and thick disk regions while obtaining 
adequate signal. This limits us to observing predom- 
inantly lower-mass galaxies, as lYoachim fc D alcantonI 
(j2006) found these are the galaxies with proportionally 
larger thick disks. 

Observations were made using the Dual Imaging Spec- 
trograph (DIS) on the ARC 3.5 meter telescope at 
Apache Point Observatory (APO) between June 2003 
and February 2006 during a total of 34 half-nights of ob- 
serving. Although the telescope aperture is much smaller 
than used for previous studies (3.5m compared to 8m), 
we were able to reach the needed signal-to-noise using a 
novel slit design. Spectroscopy of extended low surface 
brightness objects has traditionally been limited by the 
need to collect a large number of photons while main- 
taining adequate spectral resolution. However, since the 
Lick indices are defined at low spectral resolution (~ 8 
A) and cover a limited wavelength range, it is possible 
to capture more photons by using a much wider slit to 
feed a spectrograph with a higher resolution grating. We 
therefore built a custom 5"x 4' slit to allow the largest 
possible amount of light to reach the detector. On the 
blue side, we used a grating with 1200 lines/mm and 
on the red side, 830 lines/mm. The resulting spectral 
resolution was 4.9(7.1) A FWHM as measured from cali- 
bration lamps for the blue(red) chip. The blue chip gave 
wavelength coverage from 4390-5430 A sufficient to mea- 
sure the age sensitive H/3 and metallicity sensitive Mg, 
Fe5270, and Fe5335 Lick indices. The red side covered 
6140-7860 A, which allowed us to measure the strength 
of the Ha emission line. The CCDs were binned by 2 in 
the spectral direction to reduce read noise. Combining a 
large aperture slit with a high resolution grating is an un- 
conventional setup, but the increased throughput of the 
wide slit allows us to observe the low surface brightness 
regions where the thick disk dominates. 

The slit camera on DIS allowed us to accurately place 
the slit by centroiding nearby bright stars. The slit cam- 
era has a plate-scale of 0.298 "/pixel, and we were typ- 
ically able to place the slit within 1 pixel of our target 
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Fig. 1. — ij-band images of our galaxy sample with the APO 
longslit positions overlayed. All images are stretched to include 
20 < /ifl < 24.2 mag/D". 

region. Slit camera images were taken throughout long 
exposures to ensure accurate tracking. For each galaxy, 
we gathered spectra at the midplane as well as at several 
vertical scale heights above the midplane where photo- 
metric decompositions imply the majority of flux should 
be supplied by the thick disk component. Images of the 
slit placement are shown in Figure [T] and the observation 
log is listed in Table [1] 

We found that there were several prominent skylines 
near the relevant indices that would slowly fade for sev- 
eral hours after sunset. Most problematic was a skyline 
of OH at 5200 A, which often contaminated the Mg b 
index. To avoid the skylines, we took advantage of the 
APO 3.5m scheduling system which allocates observa- 
tions in half-night intervals and scheduled most of our 
observations for the second half of the night. 

3. DATA REDUCTION 

Data were processed with standard IRAF /PyRAF rou- 
tines along with several custom scripts written in IDL. 
Reduction steps for the red and blue side were identical 
unless otherwise noted. All the frames were bias cor- 
rected by subtracting the mean from the overscan region. 
Any residual bias structure was removed by subtracting 
a bias frame constructed from 5-10 bias exposures taken 
every night. For one half night (Feb 11, 2004 observa- 
tions of FGC 913), the blue chip suffered from 60 Hz noise 
which resulted in diagonal streaks across the images. 
This pattern was removed by shifting and subtracting a 
single bias frame that had similar variation in the readout 
pattern. The images were flat-fielded using spectra of a 
quartz lamp. Spectra of the twilight sky showed no need 
for an illumination correction. The spectra were wave- 
length calibrated using He, Ne, and Ar arc lamps along 



with night-sky emission lines (jOsterbrock et al.lll996f) . 
Arc spectra were taken interspersed with the science ex- 
posures throughout the night, approximately every hour. 
Stellar spectra were used to correct the spatial distortions 
(tilt) in the observations. Observations of standard stars 
and standard atmospheric extinction curves were used 
to flux calibrate the spectra. Few of the observations 
were made in photometric conditions, thus our flux cal- 
ibration is primarily used as a first order removal of the 
instrumental sensitivity profile. Because we are primar- 
ily interested in measuring equivalent widths, the exact 
flux normalization is not crucial. The sky was subtracted 
using a second order polynomial flt to regions dominated 
by the sky. The spectra were finally corrected for motion 
relative to the Local Standard of Rest, scaled to a com- 
mon flux level, spatially aligned, and combined rejecting 
cosmic ray hits. 

Our final spectra have a spatial scale of 0.42 arc- 
sec/pixel on both the red and blue chips, and a wave- 
length solution of 1.24 A/pixel in the blue and 1.68 
A/pixel in the red. Measurement of arc lamp lines 
showed a FHWM resoluti on of 4.9 A on the blue chip 
and 7.1 A on the red chip. Wor thev fc Ottavianil ()1997l ) 
report that the resolution (FWHM) of the original Lick 
indices of interest as 8.4 A. We therefore broaden our 
blue spectra with a Gaussian kernel with a = 2.9 A to 
match the Lick system resolution. 

Systematic and rotational velocities were removed by 
cross-correlating a (logarithmically-binned) stellar tem- 
plate plus Gaussian emission lines. The shifts were accu- 
rate to within 1 pixel 1.2 A). The midplane rotation 
curve was used for both the midplane and offplane spec- 
tra, as any diff erence between the two should be small at 
our resolution (jYoachim fc Dalc"anton.,2008) . 

After the 2D spectra were broadened to match the 
Lick resolution, any foreground stars were masked and 
ID spectra were extracted by summing in the spatial 
direction. For both the midplane and offplane, we ex- 
tracted spectra from the region of ±1 scale l ength (hj^), 
as measured in lYoachim fc DalcantonI (j2005[ ). 

We tested how the SNR of our spectra affects the 
accuracy of the measured Lick indices. Adding ar- 
tificial noise to template stellar spectra, we find the 
RMS error in a measured Lick index scales with S/N 
as aindex 9{SNR)^^A. Therefore, to achieve an Lick 
EW uncertainty of ±0.2A requires a SNR of 45 
per A, whic h agree s with previous determinations by 
iTrager et al.1 (|2000bD . 

Throughout our analysis, we will calculate the uncer- 
tainties in our derived ages and metallicities using the 
SNR of the extracted spectra. As an additional check, 
we have extracted ID spectra from < R/hn < 1 and 
— 1 < R/h — R < to find the variance in the Lick indices 
from one side of the galaxy to the other and note cases 
where there are large discrepancies. This procedure al- 
lows us to flag systems where systematic errors are likely 
to dominate over random uncertainties. 

4. MOVING ONTO THE LICK SYSTEM 

While we have matched th e Lick resolution reported in 
iWorthev fc Ottavianil ([1997) , this is not enough to ensure 
that we are on the Lick system. Because the Lick indices 
were originally defined from spectra that had not been 
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TABLE 1 

Edge-On Galaxy Observations 





DcitG RcingG 




Kdidpltine 




Offplcine 


Offsets 


Scale Length 


FGC 


Observed 


exposures 


total time (min) 


exposures 


total time (min) 


arcsec 


kpc^ 


arcsec 


227 


10/2004 to 10/2005 


7 


150 


20 


480 


4.2 


1.7 


10.2 


913 


02/2004 to 03/2004 


7 


80 


32 


698 


5.1 


1.5 


9.0 


1285 


02/2004 to 04/2004 


8 


105 


18 


510 


10.2 


0.9 


19.7 


1440 


03/2004 to 02/2006 


4 


80 


10 


270 


5.8 


2.0 


15.9 


1642 


06/2003 to 05/2005 


5 


100 


5 


150 


5.1 


0.9 


12.5 


1948 


06/2003 to 05/2005 


5 


95 


15 


405 


6.5 


1.2 


12.3 


2131 


06/2005 to 06/2005 


3 


45 


17 


338 


4.7 


1.0 


10.0 


2369 


10/2004 to 08/2005 


5 


105 


25 


521 


3.6 


1.0 


8.7 


2548 


10/2003 to 10/2003 


5 


75 


27 


536 


5.6 


1.5 


9.9 



^ Distances from IKarachentsev et al.l l|2000l ) 



TABLE 2 
Elliptical Galaxies Observed 



TABLE 3 
Zero point conversions 



Galaxy 


exposures 


total time (min) 


NGC 1453 


2 


25. 


NGC 1600 


2 


20. 


NGC 2778 


3 


9. 


NGC 3379 


1 


5. 


NGC 5638 


2 


10. 


NGC 6127 


2 


20. 


NGC 6702 


2 


10. 


NGC 6703 


2 


20. 


NGC 7052 


1 


5. 



.a 4 



HP 


Mgb Fe5335. 




Fe5270 




^^^^^ /^v'A/^ 








4600 4800 5000 5200 5400 6550 6600 

X (\ang) 

Fig. 2. — Examples of our blue and red spectra. From top to 
bottom: Midplane of FGC 1285, Offplanc of FGC 1285, Elliptical 
galaxy NGC 5638, Lick standard star HD114762. All the blue 
spectra have been smoothed to a resolution of 8.4 A and the Lick 
indices of interest have been labeled. 

flux-calibrated, additional corrections are needed. To do 
so, we made 144 observations of 37 unique Lick standard 
stars. Stars for which we have repeat observations have a 
mean RMS error of 0.09 A for each index . Our standard 
star E Ws are compared to values listed in lWorthev et al.l 
(fl99l . and are shown in Figure [3] The derived zero 
point corrections are listed in Table [3l In general, the 
agreement is quite good, with a typical Gaussian scat- 
ter of 0.26 A and little systematic offset. The notable 
exception is for Mg 6, for which the APO system mea- 
sures low EW for metal rich systems with Mg &> 3 A. 
We do not explicitly correct for this offset as we are ob- 
serving metal poor galaxies which have Mg b indices < 2 
A. We also comp are our standard star EWs with stars in 
common with iSchiavoni (i200Q) and find a spread of 0.2- 
0.4 A. ISchiavoni (|2006f) points out that measurements of 



index 


H/3 


Mg b 


Fe 5270 


Fe 5335 


zero point-*^ (A) 
RMS (A) 


0.13 
0.23 


0.00 
0.26 


-0.01 
0.35 


-0.05 
0.29 



^ Value subtracted from our measured indices to move 
onto the Lick system 
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0.0 

-0.5 
LO 
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2 4 

Fe5270 



-2 2 4 
Fe5335 



Fig. 3. — Comparison of Lick indices meas ured with PI S com- 
pared to the published equivalent widths in Worthev et aLl (1994). 
Dashed horizontal lines show the mean offsets while the dotted hor- 
izontal lines show the average Worthey et al. ±3(t uncertainties. 
Dashed vertical lines show the range of the line indices measured 
in our galaxy sample. Corresponding means and scatters are given 
in Table [3] 

EWs of bright stars are disturbingly inconsistent, and 
that there ca n be surpr i singly large variations between 
observations. ISchiavoni ()2006D argue that these offsets 
are caused by errors in the flat-fielding, which dominate 
the errors of the bright stars, unlike fainter galaxy spec- 
tra which are background limited. 

We also observed a sample of elliptical galaxies which 
have reported Lick measurements in the literature. The 
ellipticals were observed when light cloud cover made ob- 
serving faint regions of disks impossible, or during brief 
periods when the primary targets were not visible. We 
observed 9 galaxies in common with lTrager et all ()2000bl ) 
fisted in Table [21 Our indices are not direct ly compara- 
ble to the ones fisted in lTrager et al.l ()2000bD . as we used 
our large 5" slit which samples more flux from the outer 
regions of the galaxies, and we have not replicated the 
correction s for velocity d i spersio n and emission line fill-in 
as done in iTrager et al.l (|2000b( ). Thus, some amount of 
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Fig. 4. — Lick indices measured w ith APO compared with 
values pubhshed in ITrager et al.l Il2000bl ) . Dashed li nes show the 
median offset between our measurements and those in lTrager et al.l 
II2000H ). The agreement is good, even with our larger slit and lack 
of velocity dispersion correction. The large errors for the Fe5335 
index values are caused by the feature approaching the DIS dichroic 
cut-off. 

scatter is expected. The residuals betwee n our measure- 
ments and those in iTrager et al.l (|2000h[ ) are plotted in 
Figured! For the H/3, Mg 6, and Fe 5270 ind ices, we mea- 
sure RM S scatters of 0.16-0.31 A around the lTrager et al.l 
(|2000bf) values. There is a much larger scatter for Fe 5335 
because the feature approaches the DIS dichroic cut-off 
for galaxies with a redshift greater than 4000 km 
and becomes very low signal-to-noise. 

4.1. Emission Line Removal 

Unlike observations of gas-poor ellipticals and globular 
clusters, we must remove any emission lines that contam- 
inate the index passbands before we can derive ages and 
metallicities from the absorption features. The removal 
of emission lines is not trivial. We would like to remove 
the H/3 emission, but to do so we would need to know the 
underlying shape of the stellar continuum. This shape 
in turn depends strongly on the stellar population's age 
and metallicity, which is what we set out to measure. 
We therefore must turn to other parts of the spectrum 
to estimate the amount of H/3 emission. 

Lick indices, parti cularly H/3 , are often contaminated 
with emission lines. iGonzalea ()1993l ) popularized using 
the correction H/3emission=0.7[Oiii]A5007. This correc- 
tion was derived by fitting a stellar template to the un- 
derlying spectrum of their elliptical galaxies to isolate 
the emission feature. Unfortunately, the shape of the 
template can determine the magnitude of the emission 
correction and the ratio of f h/3 / J"foiii] is h ighly metallic- 
ity sensitive (e.g.. 'Ke wlev fc Dop ita"2002f). making this 
correction inappropr iate for low m etalli city s ystem s with 
varying stellar ages. iMoorthv fc HoltzmanI ([2006) make 
a correction by fitting a Gaussian to the emission peak 
inside the H/3 Lick index. This correction probably un- 
derestimates the H/3 correction since the Gaussian fit 
will not be sensitive to the fraction of H/3 emission that 
has filled the absorption feature. This correction is also 
very sensitive to the resolution of the measured spec- 
trum. When we try to fit a Gaussian to the emission 
feature after broadening, we find that the correction is 
0.4 A smalle r than when we m ake the correction before 
broadening. iMac Arthur! (|2006f ) takes the most extreme 
measure of masking regions where the spectrum is con- 
taminated by emission lines. 



Using the [O ill] emission line to estimate H/? is a fine 
approximation when the H/3 emission correction is small, 
as in elliptical galaxy spectra. However, our midplane 
spectra are clearly dominated by the H/? emission (Fig- 
ure [2]), and we must use a more accurate technique to re- 
move the emission. To do so, we measure the EW of the 
Ha emission line from the spectra taken simultaneously 
with the red DIS spectrograph and use this to estimate 
the EW of the H/3 emission line. 

Our procedure for removing the H/3 emission is as fol- 
lows. We first assume case B recombination with — 
2.86JHfl ([O stcrbrock 1 9891 ). This correction is u sed in 
iRampazzo et al.. (,20051 ) and iDenicolo et al.l (|2005f ) when 
measuring Lick indices in elliptical galaxies. However, 
because the Ha and H/3 emission lines were measured 
on different CCDs, we are hesitant to use the measured 
flux ratios. While both CCDs are calibrated using the 
same flux standard stars, the subsequent scaling and co- 
adding of frames was done independently and could skew 
the absolute flux calibration between the two. Instead, 
we make corrections based on equivalent widths as fol- 
lows. 

The definition of the EW is 



EW : 



1 



IX 



dX 



(1) 



where Fjx is the spectrum in the index passband and Fcx 
is the continuum spectrum calculated from the flanking 
pseudo-continuum regions. Our goal is to remove the 
contaminating emission feature to recover the equivalent- 
width of the underlying Lick absorption feature. 



EWobsorvcd = EWLick + EW, 



emission (2) 

In the case where the continuum is constant, the EW 
from the emission alone is given by 



Fcx AX 



EWemission - AA - t^l^L^Llli^^ = (3) 

i'CX i'CX 

where i^em is the integrated flux in the emission line and 
Fcx is the continuum level. We can then write the ex- 
pected EW ratio for the Balmer emission lines as 



EW(Ha)c 
EW(H/3)e 



(4) 



Fcx.na 

This in turn can be modified to account for differential 
extinction between the stars and HII regions. 



EW(Ha)c 
EW(H/3)e 



Fo.F^ 



cx,p 



lQ0AE(B-V)gasHHI3)-k(Ua) 



Fp Fcx,a 100-4S(B-l/).t„„fc(H/3)-k(Ha) 



, . (5) 

ICalzettil (|200l lists fc(H/3) - A: (Ha) = 1.163 and finds 
that E{B~V) stars = OAAE{B-V)gas- The difference in 
reddening is due to the geometrically clumpy distribution 
of HII regions compared to the smoother distribution of 
stars. Finally, we get the relation 



EW(Ha)em i^Ha F, 



EW(H/3)e 



Fcx.Ha 



CXMI3 ^q0.26E(B-V)„^ 



(6) 



For the offplane spectra, we assume there is negligible 
dust extinction, a fairly flat continuum level (-Fca,h/3 ~ 
Fcx,Ka), and case B recombination, resulting in the 
standard correction of EW(H/3)em = yEW(Ha)em with 
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y = 2.86. In the case of the midplane spectra, there is 
a younger bluer stellar population and the possibility of 
dust extinction. We adopt case B recombination, a con- 
tinuum ratio of 1.1, and E{B — V) ^ 0.1 resulting in 
y — 3.3. Our value of E{B — V) is lo wer than typical 
values for the centers of edge-on disks ([Matthews et al.l 
(|1999f) find E{B - t/) - 0.2 in an edge-on disk similar 
to our sample galaxies), however, light off the midplane 
and at larger radii experience significantly less extinction. 
For elliptical galaxies identified as dusty, iDenicolo et al.l 
(pbOSi) adopt a correction of y = 3.0. 

These adopted corrections give results consistent with 
what we expect from the broadband colors and SSP mod- 
els. Most of the corrected H/3 Lick EWs fall within 
the range expected from the model grids, and we find 
younger ages where the integrated light is blue and older 
ages where it is red. In t )4.4.11 we calculate how a 10% 
change in our adopted value of y would propagate to the 
derived ages and metallicities. This uncertainty corre- 
sponds to extinction values 0.08 < E{B — V^) < 0.42 or 
continuum ratios 1.05 < Fca,Hq/-Fca,H/3 < 1-3. 

While we consider how dust affects the emission line 
EW ratios, the Lick abso rption features th emselves are 
fairly insensitive to dust ()Mac Arthur! |2005[ ) and require 
no extra corrections even if the galaxies are dusty. 

4.2. Deriving Ages and Metallicities 

Ideally, we would fit full star formation and chemical 
evolution histories for our galaxies. However, we have 
chosen to focus on only 4 of the most prominent Lick 
indices, limiting the total number of parameters we can 
hope to fit. We therefore choose to interpret the spectra 
using SSP models. While our galaxies are clearly more 
complicated than single-burst stellar populations (given 
that we see old stellar populations along with signs of 
current star formation), the SSP models will still give 
a reasonable first-order approximation of the luminos- 
ity weighted mean ages, metallicities, and a-element en- 
hancements. 

To convert our measured indices to ages we start with 
the model grids of Thomas ct al. (2003) and interpolate 
them to a finer grid with age steps of 5t = 0.1 Gyr and 
metallicity steps of (5[Z/H]=0.01. We exclude the very 
young age grid points (age < 0.1 Gyr) as the Lick in- 
dices expected for such young ages are degenerate with 
older populations. To calculate the uncertainties, we in- 
terpolate ages and metallicities for 100 Lick index pairs 
distributed according to the equivalent width uncertain- 
ties. We also test for possible systematic errors caused 
by incorrect emission line corrections and calculate how 
a 10% error in the emission line correction for II/3 prop- 
agates to errors in age and metallicity. 

Figure O shows the Lick indices measu red for each of 
our ga laxies, placed on the model grids of iThomas et al.l 
(l200a) . assuming no a- enhancement. For each galaxy, we 
calculate the age and metallicity by interpolating onto 
grids of H/3-HMg 6, Rp+Fe 5270, and UP+Fe 5335. 

Three of the galaxies (FGC 1948, 2369, and 2548) have 
Lick indices that consistently fall far outside the model 
grids. This is not surprising, as the offplane SNR for 
FGC 1948 and 2369 were so low we should not have ex- 
pected to measure Lick indices (Table S]), and FGC 2548 
is only on the cusp of having adequate SNR. For the 
remaining six galaxies, we plot the average interpolated 



FGC1642 55kni/s 



(1 1 2 3 4 5(1 1 2 .1 ?{) 1 2 3 .1 5 

Mg b Fe 5270 Fe 5335 

FGC1642 




IZ/HI 



Fig. 5. — Lick index measurements for our observations. Top 
plots; Model grids from Thomas et al. (2003) along with our ob- 
served points. Blue points arc from midplane observations and 
red points are for offplane measurements measured from spectra 
extracted from R < Ellipses are drawn that encompass in- 

dices measured from the regions Q < R < hj{ and —hn < R < 
. Small ellipses imply the measured indices are robust, while large 
ellipses show that the spectra on either side of the galaxy are not 
particularly consistent. Bottom Plots: The interpolated ages and 
metallicities for all the above indices. Large filled points are used 
for spectra from the region —hji < R < hji that fall inside the 
model grids. The ellipses from the H/3-Mg b plane are propagated 
to the bottom plots. Throughout, blue circles are used for midplane 
observations while red triangles are used for offplane observations. 
Error bars on the points show the statistical uncertainties based 
on the SNR of the spectra. 

ages and metallicities for the midplane and offplane in 
Figure [6] and plot cumulative distributions of the age 
and metallicity differences between the thin and thick 
disk components in Figure [T] For consistency, we use 
only the metallicities measured from the Mg b and Fe 
5270 indices because several galaxies have very low SNR 
at the Fe 5335 index. 

For the 6 galaxies with high SNR, we find the thick 
disks have a median age of 7.2 Gyr and metallicity of 
[Z/II]=-0.6. The thin disks have a similar median metal- 
licity as expected for their low mass, but are uniformly 
much younger, with a median age of 4.6 Gyr. This age 
determination for the thin disk is derived after making 
very large emission line corrections and represents a flux- 
weighted average of all the stars within a radius of one 
scale length. It is certainly possible that the oldest cen- 
tral regions of the thin disk contain stars whose ages are 
similar to those of the thick disks. We return to these 
results in §4.31 below. 
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TABLE 4 

Measured Lick indices in the range —l<hii<l 



Galaxy 


SNR/Pixel 


H/3raw 


Ha 


H/3corr 




Mg b 


Fe 5270 


Fe 5335 


FGC1285 


midplane 


272.18 


-0.02±0.0 


-9.82±0.1 


2, 


.92±0, 


,0 


1.16±0.0 


1.45±0.0 


1.54±0.0 




ofTplanc 


91.21 


1.68±0.1 


-1.43±0.1 


2, 


.18±0, 


1 


1.89±0.1 


1.21±0.1 


1.27±0.1 


FGC1440 


midplane 


200.25 


-0.92±0.0 


-11.83±0.1 


2, 


.62±0, 


,0 


2.29±0.0 


1.83±0.1 


1.23±0.1 




ofTplane 


36.24 


0.59±0.2 


-4.28±0.1 


2, 


.09±0, 


,2 


2.63±0.2 


3.36±0.3 


6.75±0.5 


FGC1642 


midplane 


76.15 


0.35±0.1 


-11.28±0.1 


3, 


.73±0, 


1 


0.78±0.1 


0.68±0.1 


0.40±0.2 




offplane 


60.95 


-0.24±0.1 


-8.29±0.1 


2, 


.66±0, 


1 


1.45±0.1 


1.45±0.2 


0.84±0.3 


FGC1948 


midplane 


77.. 36 


-2.30±0.1 


-22.40±0.1 


4, 


.41±0, 


1 


0.53±0.1 


-0.16±0.1 


0.87±0.2 




oflFplane 


15.83 


0.05±0.5 


-8.99±0.1 


3, 


.19±0, 


,5 


-0.20±0.6 


-0.31±0.6 


-1.35±0.6 


FGC2131 


midplane 


98.92 


-3.82±0.1 


-29.31±0.1 


4, 


.96±0. 


.1 


0.68±0.1 


0.94±0.1 


1.82±0.1 




oflFplane 


14.71 


0.16±0.5 


-8.13±0.1 


3, 


.01±0. 


.5 


1.47±0.6 


1.51±0.7 


0.88±0.7 


FGC227 


midplane 


104.72 


-0.58±0.1 


-16.33±0.1 


4, 


.30±0. 


.1 


0.77±0.1 


1.24±0.1 


0.87±0.1 




oflFplane 


21.54 


-1.60±0.3 


-10.25±0.1 


1, 


.98±0. 


.3 


2.03±0.4 


2.45±0.4 


8.97±0.6 


FGC2369 


midplane 


51.41 


-0.70±0.1 


-14.38±0.1 


3, 


.60±0. 


.1 


1.05±0.1 


0.99±0.2 


-0.31±0.3 




ofTplane 


9.37 


1.64±0.9 


-4.84±0.1 


3 


.33±0. 


,9 


3.43±0.9 


0.19±0.9 


0.33±1.1 


FGC2548 


midplane 


81.29 


0.64±0.1 


-11.25±0.1 


4, 


■ OOiO, 


1 


-0.12±0.1 


0.99±0.1 


-0.04±0.2 




offplane 


28.38 


0.67±0.2 


-3.66±0.1 


1, 


.95±0, 


,2 


0.12±0.3 


-0.32±0.4 


-1.74±0.6 


FGC913 


midplane 


71.87 


-2.43±0.1 


-22.14±0.1 


4, 


.20±0. 


,1 


1.37±0.1 


1.30±0.2 


-0.79±0.3 




offplane 


35.48 


-0.57±0.2 


-9.45±0.1 


2, 


.73±0. 


,2 


2.86±0.2 


1.19±0.3 


-1.94±0.4 



Note. — All Lick index values are measured in angstroms. Negative values indicate emission. 
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Fig. 6. — Average ages and metallicities as measured by Lick 
indices H/3, Mg b, and Fe 5270. Blue circles show midplane (thin 
disk) observations while red triangles show offplane (thick disk) 
observations. 
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Fig. 5. — Continued. 

4.3. Radial Gradients 

The observations of the midplane of FGC 1440 were 
deep enough that we can extract radial gradients of the 
Lick indices. We extracted spectra by binning 6.3" spa- 
tially in a sliding region across the galaxy, and moving 
the observations onto the Lick system as before. This 
gives us a SNR of 200 in the central regions of the galaxy 
and 50 at 1.3 h^. The radial variation of the Lick in- 
dices along with the interpolated ages and metallicities 
are plotted in Figures [8] and H The H/3, Mg 6, and Fe 
5270 indices have radial gradients of 2.8, -1.5, and -1.3 
A/ ha respectively. These grad ients are fairly large for 
pure disk systems. iMacArth uil (2006) finds similar gra- 
dients with similar magnitudes in four late-type face-on 
systems. iMoorthy fc Holtzman (2006 ) also find strong 
radial gradients, but attribute them to the transition be- 



FlG. 7. — Cumulative distribution plots showing the differences 
between the thin disk and thick disk ages and metallicities. Pos- 
itive ages mean the thick disk is older while negative metallicity 
differences mean the thick disk is metal poor in comparison to the 
thin disk. 

tween bulge and disk dominated regions of their galaxies. 

The gradients in the Lick indices are consistent with a 
large age gradient, with the central region of FGC 1440 
having an SSP age of ^12 Gyr dropping to 2 Gyr after 
one scale length. The metallicity gradients, however, are 
slightly ambiguous. The metallicity measured from the 
Mg b index shows a mostly flat radial gradient, with a 
large drop on only one side of the galaxy. On the other 
hand the Fe 5270 index shows a decreasing metallicity 
on both sides of the galaxy. Fitting a line to the radial 
data, we find that the Mg b index has a metallicity gra- 
dient of —0.40 ± 0.1 dex/hn. If we restrict the fit to the 
region inside one scale length, the Mg b metallicity gra- 
dient is consistent with zero. The Fe 5270 index reveals 
a much steeper gradient of —0.70 ± 0.1 dex/hj^. This 
seems to imply that there is a radially changing level of 
a-enhancement throughout the galaxy, with the central 
regions being close to solar composition and the outer re- 
gions becoming more a-enhanced, thereby inflating the 
metallicity measured from Mg b. In FigurelU we plot the 
best fitting a-element enhancement and it does appear 
that the outer regions are a-enhanced compared to the 
central region which is best fit with a nearly solar chem- 
ical composition. These a-element measurements should 
be regarded with caution, as we have fit a model with 
three free parameters (age, metallicity, and a composi- 
tion) using only three measurements (II/3, Mg b, and Fe 
5270). There is also the possibility that we should adopt 
a radially varying emission line correction. This seems 
likely, as FGC 1440 hosts a dust lane which becomes less 
prominent with radius. 

While we cannot draw broad conclusions based on a 
single galaxy, it is clear we detect stronger radial gradi- 
ents than have been found in other disk systems. This 
is even more surprising given that edge-on projection ef- 
fects should act to smooth any radial population gradi- 
ents we observe. We discuss these results further in H5.5[ 

4.4. Possible Uncertainties 

4.4.1. Emission Line Correction 

Both the midplane and offplane have prominent emis- 
sion lines. For the midplane spectra, we applied II/3 cor- 
rections that had an average of 4.9 A, while the offplane 
corrections had an average of 2.2 A. If we had not take 
reddening or continuum shape into account and naively 
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Fig. 8. — Radial gradients of the Lick equivalent-widths mea- 
sured in FGC 1440. Top left shows the H/3 index after correcting 
for emission line fill-in, top right shows Mg b, lower left shows the 
Fe5270 index, and the lower right shows the uncorrected H/3 EW as 
well as the Ha EW. Dotted lines show the uncertainties calculated 
from the extracted spectra SNR. 
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Fig. 9. — The interpolated age and metallicity measured along 
the midplane of FGC 1440. The solid line shows the best fit when 
the H/3, Mg b, and Fe 5270 indices are simultaneously fit to the 
best matching IThomas et al.l I I2003I) model. Dashed and dotted 
lines show the best fit if a solar composition is assumed. There is 
a very strong age gradient present. The Mg b index shows a fairly 
flat metallicity gradient, while the Fe 5270 index shows a radially 
decreasing metallicity. The third panel shows the best fitting a- 
enhancement when we simultaneously fit all three indices. 

used just the case B recombination correction for both 
the midplane and offplane, we would have derive slightly 
younger ages for the midplane. This simple correction 
would place many of the observations near the upper 
edges of the model grids, suggesting we were under cor- 
recting. If instead we assume an even higher amount of 
dust extinction for the midplane, several of the thin and 
thick disks would then have similar SSP ages. However, 
it would take fairly extreme levels of dust extinction to 
drive all of the thin disks to have similar ages as the thick 
disks. The blue broad-band colors for these galaxies are 
inconsistent with such high levels of dust. 

Propagating an uncertainty of ±10% in the H/3 emis- 
sion correction results in a median =f0.2 dex shift in 
metallicity and ±1.9 Gyr in age for the midplane and 



=F0.1 dex and ±2.8 Gyr for the offplane. 

4.4.2. Cross Contamination 

While we have placed our longslits at regions that 
should be dominated by the thin or thick components, 
we expect some thick disk stars to be present in the 
mi dplane and vice-versa . Usi ng the photometric fits 
in lYoachim fc DalcantonI (|2006D , we find that our mid- 
plane observations typically contain ~20% thick disk 
flux, while the offplane observations have a flux contri- 
bution of ~25% from the thin disk. As can be seen in 
Figure [U we did not have large gaps between our slit 
positions. Observations made in poor seeing conditions 
might therefore experience extra cross-contamination as 
light from the midplane could be scattered into the off- 
plane slit position. This should not be a major problem 
as we avoided making offplane observations during the 
worst seeing conditions, but could slightly increase the 
a mount of expected cro ss-contamination. 

iSerra &: TrageH ()2007D study how Lick indices and their 
derived SSP ages and metallicities are affected when 
there are multiple stellar populations present. They find 
that the derived ages are very sensitive to the youngest 
stars present, while the metallicity measures are predom- 
inantly influenced by any older population. Contamina- 
tion could therefore explain the similar metallicities we 
measure for the thin and thick disk stars. However, re- 
solved stellar population studies have found small ver- 
tical metallicity gradients in low mass galaxies as well 
(|Seth et all [20051) . If enough thick disk stars contam- 
inate the midplane, our SSP derived metallicities will 
be slightly biased towards those of the older population, 
even in the young midplane. Of course, it is also possible 
that these low mass galaxies have simply not undergone 
substantial star formation episodes and thus have not 
chemically enriched the thin disks above the level of the 
thick disk. Overall, the effects of cross-contamination 
would lead us to underestimate the true metallicity dif- 
ferences between the thin and thick disks, but to overes- 
timate the flux-weighted age differences. 

4.4.3. Complex Stellar Populations 

While we are measuring SSP-equivalent ages and 
metallicities, it is fairly obvious that the midplanes of 
disk galaxies have undergone multiple epochs of star for- 
mation and are not well described by a single age and 
metallicity. Our ages and metallicities are thus best in- 
terpreted as flux-weighted averages across the extracted 
radial region, particularly for our midplane spectra where 
we have evidence for very steep radial age gradients in 
some systems. Because we are forced to bin our spec- 
tra over a large spatial region to reach adequate SNR, 
we include flux from the younger outer regions of the 
galaxy. Our midplane ages should probably be inter- 
preted as minimums, as the central region is undoubtedly 
older than the flux averaged measure we report. 

4.4.4. a-element Enhancement 

Many spectroscopic observations of elliptical galaxies 
and spiral bulges have found stellar populations that 
have systematic differences between the metallicity cal- 
culated from the Mg h index compared to the Fe in- 
dices. This systematic shift is usually interpreted as 
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being caused by a stellar population that is signifi- 
cantly enhanced with a-elements compared to the spec- 
tra that were used in building the Lick model grids. 
Such an enhancement is expected for stellar popula- 
tions that form rapidly (< 1 Gyr) an d that are primar - 
ily enriched by Type II supernovae (jMatteuccil Il994f ). 
Enhancements in g-elements are often seen in ellip- 
tical g alaxies ('Wo rthev et~all Il992l : iFisher et all Il996t 
iThoma s et al. 200^ 1, as w ell as in local MW thick disk 
stars (iBensbv et all[2005l) . 

The galaxies in our sample are all fairly low mass 
and therefore also low metallicity. In the low metal- 
licity regime, the signature of a-element enhancement 
becomes weaker in the Lick indices. Unlike massive el- 
liptical galaxies where the metallicity indicators can show 
systematic offsets of ~0.5 dex for an a-enhanced popu- 
lation, our galaxies are all sub-solar metallicity and thus 
would show little bias even if they are a-enhanced. If 
we used model SSP grids with [a/Fe]=0.3, our derived 
metallicities would change by only ~ 0.1 dex. 

We are hesitant to use our data to fit the a-element 
enhancement level. If we forge ahead and do so, we find 
considerable spread between the metallicities returned, 
but neither the thin or thick disk have systematically 
larger metallicities returned from the Mg b index, as we 
would expect if the stars were a-enhanced. However, 
given the small expected offset and lower SNR than avail- 
able for elliptical galaxy spectra, we do not consider this 
a significant result, and include it here only for complete- 
ness. 

5. DISCUSSION 

5.1. Are Low Mass Thick Disks Old? 

Measuring accurate ages for thick disk stars has been 
done in relatively few systems. In the MW, stars that are 
kinematically identified a s thick d isk stars typicall y have 
ages g reater than 8 Gyr (jFuhrm ann 1998; Bcnsb v et al.l 
l2004b[ ). HST studies of resolved stellar populations in 
nearby gala xies show the offpl ane regions are dominated 
by old stars. ISeth et al] ()2005D find that in 8 galaxies the 
offplane RGB stars have ages in the range of 2-6 Gyr. 
Similarly, iMould (2005) uses the ratio of RGB and AGB 
stars to find ages of thick disks in a sample of 4 galaxies 
to be older than 3 Gyr. Our measured thick disk SSP 
ages fall between 3.8 and 10.9 Gyr with a median age of 
7.1 Gyr, consistent with these other studies that show 
thick disks to be dominated by ancient stars. 

5.2. Are Low Mass Thick Disks Metal Poor? 

Many of our thick disks appear to be more metal rich 
than the embedded thin disks. This counter-intuitive re- 
sult is probably due to the fiux weighted nature of our 
measurement. The young, metal-poor outer regions di- 
lute the true central ages and metallicities of the thin 
disks. 

Thick disk metallicities have only been measured for a 
handful of systems. MW thick disk stars typically have 
metallicities in the range [Fe/H]~ —0.7 to -0.2, with the 
highest metallicity thick disk stars possibly reaching solar 
values (Bensby et al. 2007). One difficulty with compar- 
ing to the MW thick disk is that the observed MW prop- 
erties are of thick disk stars near the solar radius, while 
we have only been able to measure thick disk properties 



near the central regions of the galaxies. Fortunately there 
are signs that the MW thi ck disk has relativel y small age 
and metallicity gradients (|Bensbv et al.|[2005[ ). 

Constraints of thick disk metallicities have also 
been deriy e d from HST studies imaging resolved stars. 
ISeth et all (|2005h used ACS images of 6 nearby edge-on 
disks to constrain vertical gradients in the stellar popula- 
tions using the color and distribution of AGB and RGB 
stars. The older RGB stars have a systematically larger 
scale height compared to the younger AGB and main se- 
quence stars. They find little to no metallicity gradients 
in the thick disk stars in their systems, with the metallic- 
ities of the thick disks peaked around [Fe/H]~ —1. This 
is slightly more metal poor than the measurements we 
have for our thick disks. Howe ver, the offset is u nlikely 
to be significant, given that the ISeth et aH (|2005f ) study 
is able to study a cleaner sample of thick disk stars by 
reaching higher vertical heights which reducing the con- 
tami nation of thin dis k stars. Like the data presented 
here, Seth et all ()2005l ) only studies lower mass systems, 
limiting the amplitude of any possible metallicity gradi- 
ent due to the low metal hcity o f the m idplane. Using 
similar HST observations, iMouldl (|2005( l finds that thick 
disk stars in 4 edge-on galaxies have [Fe/H] between -1.0 
and -0.78, again very similar to the metallicities we find. 

In Figure [TOl we comp are our t hin disk metallicities to 
the low mass sai nple in iLee et a l. (2006) and the large 
SDSS sample of iTremonti et al.l (,2004) . We also com- 
pare our thic k disk values to th i ck disk and halo s ample s 
presented i n 'Seth et all (|2005'): 'Mouhcine et al.' (20051); 
iTikhonov e t al. (200'^andFRcddv et al. (2006). Incases 
where the authors presented values of log(0/H)+12, we 
converted to [Fe/H] assuming log(O/H)0-f 12=8.69 and 
[Fe/O]=0. We also plot the midplane nebular abun- 
dan ces for our g alaxies calculated from the S2N2 calibra- 
tor (jViironen et al. 2007). With the exception of a few 
outliers, our measured metallicities are consistent with 
metallicities measured in simila r systems. 

Recently, llvezic et"al] (|2008f ) have questioned if the 
Milky Way thick disk is a unique component or simply an 
extension of the thin disk that has non-Gaussian metal- 
licity and velocity distributions (|Norrislll987[ ). They cite 
a lack of correlation between velocity and metallicity in 
large SDSS samples as a major problem for traditional 
disk decompositions and rule out a "traditional" two- 
disk model at the 8a level. Unfortunately, the model 
llvezic et al.l (|2008l ) rule out is not applicable to their 
observations. In particular, they model 1,1 42 stars ob- 
served in the region 1.0 < z/kpc < 1.2. llvezic et al.l 
()2008l ) correctly assume the observations will contain a 
similar number of thin and thick disk stars, however they 
do not use a realistic thin disk component. The thin disk 
stars are modeled as having an asymmetric drift of 9 
km s~^and [Fe/H]=-0.50 with a spread of 0.04 dex (their 
Figure 16). These parameters would be appropriate for 
modeling nearby thin disk stars, but they are observing 
stars 3-4 thin disk scale heights above the Galactic plane. 
At this large height, only the kinematically hottest thin 
disk stars will be present in the sample, and one should 
expect a much larger velocity lag than observed in local 
cooler thin disk stars. iHolmberg et al.l ()2007l ) show that 
the hottest local disk stars are the oldest and that the 
age-metallicity relation (AMR) for local stars has large 
intrinsic dispersion at large ages {a ^ 0.2 dex). This 
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Fig. 10. — Comparison of our thin and thick disk metallic- 
ities to similar studies. The panel on the left shows our thin 
disk metallicitios along with the low mass gal axies in |Lee et al.l 
J2OO6) and the fits to larger samples in Trcmon ti et al.l (120041 ) and 
lYoachim & Dalcanton 1 20071 ). On the right, we compare our thick 
disk metallicitios to the systems me a sured in [Sctli et al. (2005), 
[Mouhcine et al.i I i2005l) . IReddv et all (1200^) . and iTikhonov et al.i 
II2005I) . 

large metallicity dispersion is probably a result of ra- 
dial migration o f the older high velocity dispersion stars 
(|Havwoodll2008l) . 

If the offplane thin disk stars are assumed to be simi- 
lar to older local thin disk stars (i.e., kinematically hot 
with a broad metallicity distribution and larger asym- 
metric drift), their properties should well match SDSS 
observations of regions at large scale height. Specifi- 
cally, there should be little correlation between kinemat- 
ics and metallicity despite the mixture of thin and thick 
disk stars. This revised model is consistent with the ob- 
servations o f thick disks presented here and elsewhere 
(jSeth et al.l[200 5) that find stars at large scale heights 
are a significantly older population than those found near 
the midplane. 

5.3. Are Thick Disks a-enhanced? 

The level of a-enhancement can be a major clue to the 
formation process of a stellar population. If a-enhanced, 
it is a sign that a stellar population has been enriched 
mostly over a short time period by Type II SNe, whereas 
stars with solar composition formed over an extended pe- 
riod and have been enriched by both Type II and Type 
la SNe. Numerous papers have found that MW thick 
disk stars are enhanced in a-elements co mpared to thin 
disk s t ars at s i milar t otal metallicities (e.g.jBensbv et al.l 
20pl d^Ial, I2OO5I; iFeltzing et all [20o| I Reddv et all 
2006; Tautvaisicnc ct al. 2001: Mashonkina et al l l200a 
Prochaska et al.,,200Q : ,Fuhrmann..l99&. ,2004i) . Unfortu- 
nately, we do not have the SNR to definitively say if our 
thick disks are a-enhanced, due to the overall low metal- 
licities of our target galaxies. 

5.4. Are Thick Disks "Normal" Stellar Populations? 

Using a sample of >1000 S PSS images of edge-on disk 
galaxies. IZibetti et al.l (|2004f ) examine the faint halo that 
appears when the images are stacked. This extended 
halo has anomalous colors, requiring stars t hat are either 
metal ric h or have a bottom heavy IMF (jZibetti et al.l 
(2004; Zack risson et al.l l 2006f ). We find no such anomalies 
with the thick disks we observe spectroscopically, as most 
of them fall on stellar synthesis model grids using stan- 
dard IMFs and metallicities. The few galaxies where we 



measure Lick indices that are inconsistent with the SSP 
models can easily be explained as spurious measurements 
caused by low SNR, and do not require exotic stellar pop- 
ulations. 

5.5. Radial Color Gradients in the Thin Disks 

iBell fc de Jond (|2000D observed broadband colors for 
a large sample of galaxies and found that the radial gra- 
dients were predominantly caused by age gradients in 
the stellar populations. By averaging their sample to- 
gether. They found a metallicity gradient, but broad- 
band sensitivity to dust makes this a measurement in 
in dividual galaxies. Their work has been expanded on 
bv lMacArthur et all (|2004D . who found both mctaUicity 
and age gradients are stronger in the inner regions of 
galaxies, and that galaxies with strong age gradients had 
smaller metallicity gradients. The Lick indices ability to 
lift the age-metallicity degeneracy, and relative insensi- 
tivity to dust makes it much easier to quantify how much 
of the radial color gradients in disk galaxies are due to 
age versus metallicity changes. 

In one of the only study that has explicitly targeted 
Lick indices in th e disk-dominated regions of galaxies, 
iMacArthuii ()2006( ) detect age gradients in only 2 of their 
8 galaxies, and find strong negat ive meta llicity gradients 
in 4. The age gradients in , Mac Arthur! j2006.) are also 
rather small (-0.5 and -1.3 Gyr//i/j). Our measurements 
of FGC 1440 show a much steeper age gradient, with the 
S SP age dropp i ng by 9 Gyr over one scale length. 

iGanda et all (|2007f) use the SAURON integral field 
unit spectrograph to measure Lick index strengths across 
the face of 18 late-type disk galaxies. They find in general 
up increasing with radius and the metal sensitive lines 
decreasing with radius. The SAURON observations show 
a very wide range of galaxy-to-galaxy radial behavior, the 
most extreme are consistent with the strong gradients we 
find in FGC 1440. 

Our finding that the thin disk of FGC 1440 might 
have stellar populations of near solar composition and 
be a-enhanced at larger radii is puzzling. Looking at 
the t ransition between bulge and disk dominated re- 
gions, [MQorth^5!Soitimaa([2006|) finds that the central 
bulges are either solar-composition or a-enhanced, with 
little to no a-enhancement in the disks. With the pres- 
ence of old thick disks in all of the galaxies, we would 
expect the galaxies to have experienced plenty of chemi- 
cal evolution and enrichment from SNe la. Instead, the 
outer regions of FGC 1440 are a-enhanced, suggesting 
that the central region of the galaxy has undergone ex- 
tended chemical enrichment, while the outer regions have 
not, despite being surrounded by old thick disk stars. 

This could be a sign that the thick disk stars in FGC 
1440 have been recently accreted, and thus have not con- 
tributed to the chemical enrichment the galaxy. An- 
other possibility is that the central region of the galaxy is 
the only place where the gravitational potential is deep 
enough to retain SN ejecta, and the outer disk has histori- 
cally suffered from SN blow-out and failed to retain metal 
enriched gas. Another possibility is that the luminosity- 
weighted metallicity in the outer disk is dominated by 
enrichment from the latest burst of star-formation which 
has a-enhanced the region. 

6. CONCLUSIONS 



Lick Indicies of Thick Disks 
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We have spectroscopically confirmed that the thick 
disks observed in edge-on late type galaxies are old, 
metal-poor stellar populations, analogous to the thick 
disk stars seen in the MW and nearby edge-on systems. 
This is the first time ages and metallicities of thick disks 
have been measured in unresolved stellar populations. 
Because all of our targets are fairly low mass, we are 
unable to detect any significant differences between thin 
disk and thick disk metallicities. After correcting for 
emission line contamination, the thin disks in our sam- 
ple are found to be quite young, with strong radial age 
gradients. 

We fail to detect any significant trend for thick disk 
stars to be enhanced in a-elements compared to their 
thin disks which is a defining characteristic of the MW 
thick disk. Our failure to observe a-enhancements is 
most likely a result of our sample being dominated by 
low-mass and therefore low metallicity galaxies, for which 
solar and a-enhanced models are similar. 
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